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Abstract 

We study the cosmological consequences of generic scalar fields like moduli 
which decay only through gravitationally suppressed interactions. We con- 
sider a new production mechanism of gravitinos from moduli decay, which 
might be more effective than previously known mechanisms, and calculate 
the final gravitino-to-entropy ratio to compare with the constraints imposed 
by successful big bang nucleosynthesis (BBN) etc., taking possible hadronic 
decays of gravitinos into account. We find the modulus mass smaller than 
~ 10^ TeV is excluded. On the other hand, inflation models with high re- 
heating temperatures jnf ~ 10^^ GeV can be compatible with BBN thanks 
to the late-time entropy production from the moduli decay if model parame- 
ters are appropriately chosen. 
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Typeset using REVT^ 



I. INTRODUCTION 



One of the consequences of local supersymmetry or supergravity is the existence of grav- 
itinos, the superpartner of gravitons, whose natural mass scale is the weak scale 0(0.1 — 10) 
TeV. If they are not the hghtest supersymmetric particle (LSP), they decay into LSP and 
other high energy particles such as photons, neutrinos, quark-antiquark pairs or gluons af- 
ter/during the big-bang nucleosynthesis (BBN) epoch. Such decay products may change 
the light element abundances by changing the baryon-to-entropy, neutron-to-proton ratios, 
or destroying the produced elements, which may result in a significant discrepancy between 
yields and observations. This is called the "gravitino problem". A lot of authors have 
studied a variety of constraints on the decaying gravitinos from BBN [1-23] ^. 

In the standard inflationary cosmology, gravitinos are produced by the scattering pro- 
cesses of thermal particles at the reheating epoch just after inflation. The yield parameter 
?/3/2 of gravitinos, which is the ratio of their number density to entropy density, 1/3/2 = 
is approximately given as a function of the reheating temperature after inflation, TR^mf, as 
[15,27,28] 2 

.3«.-1.5xl0-(^). (1) 

Hence it is customary to express the constraint on their abundance imposed by BBN as 
that on the reheating temperature Ts^nf- For example, if gravitino mass is equal to 0.1 TeV 
we find an upper bound T^jinf ~ lO^GcV, which imposes a constraint on model building of 
inflation. For example, hybrid inflation [29] is difficult to reconcile with this low reheating 
temperature, in which the inflaton is typically coupled to gauge fields and decays rapidly. 
So it is preferable to have another inflation after hybrid inflation as suggested by recent 
observational data [30] . 

Wc would like to point out, however, that these constraints on T^^nf have all been 
obtained under the implicit assumption that y3/2,inf remains constant until their lifetime, 

where we assume that the gravitino decays into a massless gauge boson and a gaugino, A^g 
is the number of the generators of the gauge group, and m3/2 is the gravitino mass. In this 
paper we argue that both the denominator and the numerator of 1/3/2 = are subject 

to change between the reheating epoch after inflation and their decay time, apart from the 
dilution due to cosmic expansion which does not change the ratio 1/3/2 itself. We then derive 
more appropriate constraints imposed by BBN. 

In supergravity or superstring theories there appear a number of long-lived scalar flclds 
which decay only through gravitational interactions, such as moduli, dilatons, or Polonyi 



"^Cosmological constraints on stable gravitinos are studied in Refs. [24-26] and references therein. 

^The error of this approximation formula is within ~ 5% for Tr = 10^ - 10''^^ GeV, and ~ 25% 
for Tr = 10^ - 10^9 GeV [27]. 
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field, which are referred to collectively as the modulus hereafter. The modulus (f) starts 
coherent field oscillation to dissipate its energy density as the Hubble parameter H becomes 
smaller than its mass m^. Because its dissipation rate is smaller than that of radiation, the 
universe turns to be matter dominated well before the lifetime of (f), when it decays producing 
not only huge amount of entropy, causing what is called the moduli problem [31-33], but also 
undesirable particles for cosmology. So far a number of groups have studied cosmological 
constraints on the modulus decays depending on the properties of each decay product, e.g., 
for LSPs not to close the universe [34-36], and for radiations to complete thermalization [37]. 

The effects of decaying moduh on the gravitino problem are twofold. One is that the 
entropy produced by their decay dilutes primordial gravitino abundance, which is a good 
news to relax constraints on inflation model building. The other is that these unwanted 
particles may also be produced directly by the decay of moduli (f). Including these two 
effects, it is not apparent for us to find the allowed region of model parameters to avoid 
both the gravitino problem and moduli problem. So far no one has considered this type 
of scenario with decaying moduli. Therefore, we comprehensively study the effects in this 
work. 

The constraint by gravitinos produced directly by the decay of moduli was first studied 
by Hashimoto et al. [38], who considered the case that the modulus decays into two grav- 
itinos. In the present paper, however, we point out a more efficient mechanism of direct 
gravitino production from moduli decay, which is the mode that (p decays into its superpart- 
ner modulino, 0, and a gravitino, V'3/2, namely, + '03/2- We incorporate this decay 
mode and consider cosmological constraints on the masses of gravitinos and moduli from 
BBN and other observations. 

The fate of the decaying moduli depends on its mass m^. In [32,33], it was presumed 
that moduli fields acquire masses through supersymmetry breaking and then their masses 
are comparable to masses of superparticles. In this case the modulus field is long-lived 
and the reheating at the modulus decay takes place with the temperature much below 
IMeV, spoiling the success of the BBN. (Sec, however, [34]). Recently it was realized that a 
mechanism to stabilize the moduli fields is operative in the compactification with non-zero 
NS and RR fiuxes in certain string theories [39]. This makes most of the moduli fields 
very massive, typically around the string scale. Still there are some moduli which are not 
stabilized. In particular, the modulus field which is responsible to determine the size of 
the compactification is not stabilized in the flux compactiflcation. With the ignorance of 
possible mechanisms on mass generation, we take the mass of the modulus field as a free 
parameter in this paper. 



II. MODULUS DECAY INTO GRAVITINOS 



The relevant terms in the supergravity Lagrangian which describe the decay mode of our 
concern are given by 



1 



G 



(3) 



This is allowed by gauge invar iancc. Furthermore it has the same structure as a source of 
soft supersymmetry (SUSY) breaking terms of squarks and sleptons in gravity mediation 
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of SUSY breaking (see e.g. [40]). Here a^^^^ is the dimensionless coupling constant, Z — 
z + \/2dil>z + OOFz + ■ ■ ■ is the chiral superfield which breaks SUSY by acquiring an F term. 
Therefore, it includes goldstino ipz or the longitudinal mode of gravitino, with a complex 
scalar field z and an auxiliary field while 9 and 9 are two-component Grassmann variable. 
$ = (f)+\/29^+99F^+i^/29dn^9a'^9+- ■ ■ is the chiral superfield which includes the modulus 
0, and the modulino cp with an auxiliary field F^. Then we obtain 



jC — a. 



^4, 



(4) 



where we have incorporated the SUSY breaking effect, {F^ = \/3ms/2MG with the reduced 
Planck mass Mq (= 2.44 x 10^^ GeV) without generating a cosmological constant. We also 

used the field equations of motion, for the modulino i(T^9^0 — m^cj) — with its mass m^, 

and for the modulus F^ + m^(f) = 0. 

Then, the mass matrix of is given by 



mi 



(5) 



Therefore, we can diagonalize it and get two eigen values m</,_,_, which are 



± 



mi 



(6) 



We find that the decay ^ + V'3/2 is kinematically allowed if a. 



> 



(2 + ms/2/m^ 



In this paper, we consider the case m^_^ — -C or ^ y '-"^iIjzH 



> \/3a^^00 so that the 

field oscillation of 0+ continues as long as that of 0_. Hereafter for simplicity we assume 
that only 0+ is present and induces a coherent oscillation with mass. 



~ + -Y^^^zU "^3/2- (7) 
As is seen in (4), the relevant term of the Lagrangian that describes the decay reads 



mi 



G 



(8) 



Then the decay rate of the modulus into the gravitino ipz and the modulino is given by 



Stt 



m^ 



( X 2-1 3/2 r 



m3/2 



21 1/2 



(9) 



Because we are assuming m<^ ^ V^o^,^^^ ?7t.3/2, we obtain the following approximate formula 
for Eq. (9), 



4 



2 




(10) 

where the factor ^4 is given by 

"^^^ -1) (11) 

Note that this decay rate is very sensitive to the couphng constant a^^^^- Therefore, for 

relatively larger a^^^^, the decay width of the mode — > ■03/2 + can become larger than 
that of the mode (p 2ii>3/2 which was studied in Ref. [38] ^. 

On the other hand, the decay width of modulus into radiation is represented by 

■) 

r(0 — >• radiations) = ^ 2 1 (1"^) 

where N depends on the number of the final states. For example, if all of the particle 
contents in the minimal supersymmetric standard model (MSSM) appear in the final states, 
we approximately obtain N = (9(10) [36]. Because we are interested in the parameter space 
which satisfies ^ v^o,^^ ^^7723/2, the above decay width into radiation is the dominant 
decay mode, which we identify with the total width, Ttot hereafter. Thus the branching ratio 
to gravitino production reads 



III. GRAVITINO ABUNDANCE AFTER MODULUS DECAY 

Now we consider cosmological evolution of the modulus (p. By adding an appropriate 
constant we redefine so that it has a global minimum at = and assume that the mass 
term dominates its potential energy density for simplicity. We also assume that its initial 
amplitude, (pi, is of order of Mq or smaller. The modulus remains there until the Hubble 
parameter H decreases to m^, when it starts coherent oscillation around the origin. As the 
field oscillation red-shifts less rapidly than radiation, the universe will be dominated by 4> 
at the time. 



^The decay width of the mode cf) 2^^3/2 is proportional to the square of the relevant coupling 
constant as opposed to oc here. On the other hand, if the coupling constant is smaller than 

unity or the decay mode (p — ^ "03/2 + f?^ is kinematically forbidden, the decay mode (p 2^3/2 is 
more effective. 
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After that time, the expansion law of the cosmic scale factor a(t) is the same as that in 
matter dominated regime, a(t) oc t^/'^, until (f) decays at td ~ F^t, when the ratio of energy 
density of 4> to that of radiation reads 



PR{td) \teq) 



9N \ J \Mg^ 



(16) 



Thus the entropy increase factor, A, reads 



Sbefore(^d) \PR(td) J ^ 9N 171^ , 



(17) 



where Safter(^d) ( =^bcforc(i^d) ) IS the entropy density just after (before) the modulus decay. In 
this paper we only consider the case that A 3> 1. 

Using an approximation that the modulus energy density is fully converted to radiation 
when Ftot = H{Tji), the reheating temperature after the modulus decay is found to be 

Here denotes total effective numbers of relativistic degrees of freedom. We find = 10.75 
just before the onset of BBN and g^ = 228.75 if all the particle contents of the minimal 
supersymmetric standard model are massless and in thermal equilibrium. We note that the 
reheating temperature should satisfy 

Tr > 1.2 MeV, (at 95% C.L.), (19) 

so that the neutrino background can complete thermalization to warrant successful 
BBN [37] ^.This means that should satisfy^ 



m^>117V-V3(^_|^y^'TeV. (20) 



^In Ref. [37], however, Tr is defined by Ftot = 3-ff(TR). Using their definition, the lower bound 
of Tr turns into 0.7 MeV, which they reported. On the other hand, note that they also discussed 
constraints for emitted hadrons by decaying modulus not to influence on the neutron to proton 
ratio before/during the BBN epoch. Then, the lower bound is pushed to the severer one {Tr > 5.2 
MeV). In this paper, however, we do not go into such specifics. Here we adopt the conservative 
one (Tr > 1.2 MeV). 

^This constraint may be evaded if the Universe underwent late-time inflation such as the ther- 
mal inflation to dilute the energy density of the modulus field [41]. In this case, the reheating 
temperature after the thermal inflation should satisfy this limit. 
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Now there are three sources of gravitinos after the modulus decay. One is the primordial 
gravitinos which were produced just after inflation and diluted by the entropy from moduli. 
Using (1), its abundance is given by 

. . y3/2,inf (tj) 1 , in-21 Arl/2 f '^4> \ f ^l^i \ ^ f J^RM_\ x 

y^/^MiU) - - ^-^ ^ N [WfeVjyM^) llO^^j' ^^^^ 

where ?/3/2,inf (^i) refers to the primordial value after inflation and we have assumed no sig- 
nificant entropy production took place between reheating after inflation and modulus decay. 

On the other hand, the yield parameter due to direct production of gravitinos from (f) is 
given by 



= 1.1 X 10 



-16 



3A4 






SttN 










( 9* 


3 J 


V102 



UTevy uoexevy ■ ^ ^ 



Here we have used the following approximate relations; s = (27r^/45)5'^,T^, and pji — 
(7r^/30)(7*T^ with the energy density of radiation at the reheating time. 

Finally, gravitinos are also produced by the scattering process in the thermal bath at the 
reheating due to modulus decay, whose contribution to the yield parameter reads 



y./2Ata) ^ 1.5 X 10-^^ (yo^ 



= 1.8 X 10 



-18; 



Comparing these three equations with each other we can find which production mech- 
anism is dominant for each combination of model parameters. We first calculate which of 
the late time production mechanism is more efficient. From (22) and (23) we find the direct 
production from decaying moduli is more efficient than the thermal scattering processes if 

< 4 X 10^ (^)''' N-y^ (^)''' TeV ^ m^.. (24) 

Then comparing y3/2,th(^d) with y3/2,inf (^d) for the case > m^ci and yz/2,4>{td) with 
y3/2,inf(^d) for the opposite case < m^cr, we find that primordial gravitinos can dom- 
inate over the late-time counterparts only if the inequality 

is satisfied, where the factor J is determined by the combination of efficiency of primordial 
production of gravitinos and deficiency of dilution due to entropy production from (p. 

First suppose that the inequality (25) is satisfied. Then we find that the direct production 
from (h is dominant for 



%.^A'^' (9*_\-'"' (!!^" (J_\-"\ 



m,<2xlOM^) [j^) TeV, (26) 
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while primordial one contributes the most for 

\ 3 ) VlOV VlTeV/' VlOV 



<m,<6xlO'(|,)"^(^)\eV. 



(27) 



For the case 



the thermal scattering in the plasma produced by modulus decay is the most important. 
The situation with J < is much simpler. We find that y^/2,(i>iid) is dominant if 
< m^cr and that yz/2,thiid) is the largest for > m^cr- 

IV. COSMOLOGICAL CONSTRAINTS ON DECAYING GRAVITINOS 

Having fully analyzed the yield parameter for all the possible cases in our model, we are 
now in a position to impose cosmological constraints on our model parameters. First we 
summarize the constraints we use. 

If the gravitino further decays into lighter particles such as photons, neutrinos, gluons, 
or quark-antiquark pairs, they induce electromagnetic or hadronic showers which would 
influence on BBN because they may change the light element abundances by destroying 
them or changing the neutron to proton ratio. So far a number of authors have studied a 
variety of constraints on decaying gravitinos from observations [1-23] . Here we consider the 
following two typical examples. One is the case that the gravitino decays only into a photon 
and an LSP neutralino x with Nq — 1. Even in this case, there also exists a hadronic decay 
mode into quark-antiquark pairs with the hadronic branching ratio Bh ~ a/{47r) ~ 10~^. 
The other is that the gravitino decays only into a gluon g and a gluino g with Nq = 8. 
Then, the hadronic branching ratio is unity Bh = 1. In both cases, the electromagnetic 
branching ratio is practically unity because almost all energies of their decay products are 
immediately turned into high energy radiations that participate in the photodissociation 
of light elements. Including more hadronic decay modes in addition to the radiative decay 
mode, it is known that the resultant constraints get more stringent [11-13,21,23]. 

Compared with observational light element abundances, wc get the upper bound on 1/3/2 
as a function of the gravitino mass. It is approximately expressed by [21,23] 

3 X 10-1^ - 6 X 10"^^ (3 X 10-18 - 2 X IQ-^^) for 7713/2 = 0.1 - ITeV, 
2 X 10-15 _ 2 X 10-^^ (3 X 10-18 - 2 X 10-13) for 7713/2 = ITeV - lOTeV, (29) 
2 X 10-13- ~ lO-'^ (1 X 10-13 - 2 X 10-13) for 7773/2 = lOTeV - 30TeV, 

for Bh — 10-3 (Bh = 1)- Of course, if we do not consider the decaying moduh, (29) is 
transformed into upper bounds on the reheating temperature after the primordial inflation 
Tn^mi as follows. 

2 X 10^ - 4 X 10^ (2 X 10^ - 1 X 10^) for 7773/2 = 0.1 - ITeV, 
7ii,inf < <( 1 X 10^ - 1 X 10^ (2 X 10^ - 2 X 10^) for 7773/2 = ITcV - lOTeV, (30) 
1 X 10^ 1013 (1 X 10^ - 2 X 10^) for 7773/2 = lOTeV - 30TeV, 



2/3/2 ^ 
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for Bh — 10~^ {Bh = !)• Satisfying these upper bounds is just one solution to avoid the 
gravitino problem without any late-time entropy production by decaying particles such as 
moduli. For more detail, see the results in Ref. [21,23] ^. 

In the present scenario, an LSP is produced by each gravitino decay. Therefore we should 
consider the constraints on relic density of LSPs. We can relate their number density riLSP 
with that of gravitinos, 

y3/2 = = , 31) 

s mLSP s 

where mLSP and Qlsp are the mass and the density parameter of LSPs, respectively, and pcr 
is the critical density in the universe.^ Prom observation of CMB anisotropics, the WMAP 
collaboration reported that JIlsp < 0.35 at 95 % C.L. [42]. Here, p^/s is equal to the 

present value Pcro/so, where Pcro = 4.2 x 10-^^GeV^(V0.72)2 with h = 0.72 ± 0.05 [42,43], 
and So = 2.2 x 10~^^GeV^(To/2.725K)^ with the present value of the photon temperature 
To = 2.725 ±0.002K [44]. Normalizing at mLSP = 100 GeV, we get the upper bound on 1/3/2 
from Eq. (31), 

which is a conservative bound that turns out to be exact in the case thermal relic density is 
negligibly small.® Here we have assumed that the abundance of these neutralinos produced 
by such decaying gravitinos does not change though possible annihilation process, which can 
be shown to be ineffective as follows. The annihilation cross section of the neutralino a^nn 
is a complicated function of the superparticle mass spectrum, but is generically bounded as 



ml 



(33) 



with the coupling ~ 10 ^. Then, the ratio of the annihilation rate Pann = ^T-LSpCann, to 
the Hubble parameter H ~ T'^/Mq is bounded as 



^The radiative decay of gravitinos also influence on the shape of Planck distribution of cosmic 
microwave background (CMB) radiation through the fi or y distortion. However, this type of limit 
is weaker than that of the photodissociation process [14,18]. 

our scenario, the modulus decay also produces a modulino, followed by the decay into a LSP. 
Thus the total amount of the LSPs could be up to twice as large as the right-hand-side of Eq. (31). 
We do not include this effect in our analysis, which does not affect our conclusions. 

^In fact, as the modulus mass increases, the reheating temperature exceeds the freeze-out tem- 
perature above which the thermal production/annihilation of the LSPs is effective. In this case 
a significant contribution to J^lsp may come from the thermal production soon after the moduli 
decay. This contribution is very sensitive to the SUSY mass spectrum and, for simplicity, we do 
not include it in this paper. 
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We find Fann/ H < 1 at the temperature of gravitino decay. 

Tlie LSPs may also be produced by the direct moduhis decay into the superparticles in 
the MSSM. This issue was investigated in Ref. [36]. It turned out that the result depends 
on the couplings of the modulus field to the MSSM fields. Consider, for instance, the case 
where the modulus field couples to gauge multiplets through gauge kinetic function. Then 
the modulus decay into a pair of gauginos receives chirality suppression, and the branching 
ratio is suppressed by {rrig/m^Y, with rrig being the gaugino mass [36]. A rough estimate of 
the number density of the LSPs produced by this decay gives 

^^Br(^^m^~(^y(^)"\ (35) 

s 777,0 \m^J \MgJ 

which is negligibly small for the range of the mass parameters of our concern. Inclusion of the 
moduli coupling to quark/lepton chiral multiplets in Kahler potential does not change the 
result, because the latter decay is further suppressed by the fourth power of squark/slepton 
masses. In this paper, we assume that this is the case, and we do not consider the direct 
production of the LSPs by the modulus decay. From the above discussions about LSP, one 
can see that our treatment gives conservative limits. 

We now depict our constraints on 7773/2 and 777^. First in order to clarify the influences 
on cosmology from the gravitino produced only by the modulus decay, i.e., from y3/2,(j,{td) 
and y3/2,th(^(i)j for the moment we assume that ?/3/2,mf is negligibly small. This condition is 
represented by J < J^r- If the reheating temperature after inflation T^inf is not so high, 
i.e., ^ 10^^ GeV, this situation is realized with 0^ ~ Mq. Then yz/2,\ni is entirely diluted 
by the late-time entropy production caused by the moduli decay. Using the observational 
upper bounds in (29) and (32), we can get constraints on modulus mass 777^ as a function of 
gravitino mass 7773/2. From the expressions of 1/3/2 (^d) in Eq. (22) and Eq. (23), we can get 
the lower and upper bounds on m^, respectively. 

In Fig. 1, we plot the resultant constraints on 777^ as a function of 7773/2 for the case 
l/3/2,inf = and the hadronic branching ratio Bh — 10"^. Here we have taken a^^^^ ~ 3 
and N — 1 in Eq. (13). We find that the reheating requirement (19) gives a milder lower 
bound on m^. Prom Fig. 1, we see that the modulus mass of the weak scale is excluded for 
gravitino mass of 7773/2 = 0.1 - 100 TeV. In addition, it is interesting that we can obtain 
the upper bound on 777^ by this type of cosmological arguments. Finally we comment on the 
dependence of our constraints in Fig. 1 on N. Since appears only in the form N^l'^m^l'^ 
in all of the relevant expressions in Eqs. (18), (22) and (23), the constraints for N other 
than iV = 1 can easily be read off by replacing 777^ in the vertical axis by N^^^m^, in Fig. 1. 

Next we discuss the more general situation where 7/3/2,inf can also be important with 
J ^ Jcr. We plot the constraints on 777^ as a function of 7773/2 in Fig. 2 in the case Bh = 10~^. 
The extent of oblique lines coincides exactly with those excluded from the constraints by 
BBN and LSP in Fig. 1, which correspond to the hmiting cases that TR^mf is sufficiently low 
such as T^^^ = — 10^^ GeV, namely, J < J^. For higher reheating temperatures e.g., Tjii^f 
= 10^'^ - 10^^ GeV, larger parameter regions are additionally excluded. The shadowed region 
corresponds to the excluded region for Ts^inf = 10^^ GeV with 4>i = Mq or for J = 10^. In 
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a similar fashion, the constraints for the case Bh — 1 are depicted in Fig. 3 and Fig. 4 for 
J < JcT and J ^ Jcr, respectively. From these figures, we see that larger regions are excluded 
for Bh = 1. 

It is important to note that although the excluded region becomes broader for higher 
reheating temperature, we still have a fairly large allowed region in our parameter space 
even for the highest possible reheating temperature T^^inf = 10^^ GeV thanks to the dilution 
of primordial gravitinos by the entropy production associated with modulus decay. Thus 
the previous upper bound on Tr mf in (30) can easily be relaxed if we consider the decaying 
moduli. In Table 1 (Table 2) we show the allowed values of m<^ for various 777.3/2 and TR,inf 
for the case Bh = 10"=^ {Bh = 1). 

V. CONCLUSION 

Wc have studied the effects of decaying modulus oscillation on the cosmological gravitino 
problem. We have considered a new direct production mechanism of gravitinos from modulus 
decay, namely, a decay mode of modulus into a gravitino and a modulino. The width of this 
decay mode can be larger than the other mode into two gravitinos which has been studied 
in [38], if the coupling constant is of the same order of magnitude with n,,^,^^^ > 1. 

Comparing our yield of gravitinos with the constraints imposed by BBN and the relic 
LSPs, which are decay products of gravitinos, we have obtained a constraint on the masses of 
gravitinos and modulus. As a result we have found that due to the above-mentioned direct 
production of gravitinos from decaying modulus, the modulus mass with 777^ < lO^TeV is 
excluded, even when the branching ratio into hadrons is minimal. 

On the other hand, wc have also found that wide range of 7/73/2 and 777,^ are still allowed 
even if the reheating temperature after inflation is as high as TR,i„f = W GeV and the 
effects on the hadronic decay of the gravitinos are taken into account, thanks to the dilution 
of primordial gravitinos due to the entropy production associated with modulus decay. 

Thus in order to study cosmological consequences of gravitinos, it is important to analyze 
not only their abundance right after inflation but also their subsequent dilution due to late- 
time entropy production, as well as late-time production from scalar condensates with only 
gravitationally suppressed interactions including a dilaton and a Polonyi field. 
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FIGURES 




FIG. 1. Constraints on modulus mass m,^ as a function of gravitino mass w.3/2 in case that 
the hadronic branching ratio is Bh = 10~^. The sohd hues denote the hmits which come from 
the photodissociation and hadro-dissociation effects on hght elements by the decaying gravitino 
(BBN). The dotted lines represent constraints on the energy density for LSPs, which are produced 
by decaying gravitinos, to satisfy J^lsp < 0.35 for m-LSP = 100 GeV. The dashed line denotes 
the lower bound on m,^ arising from the lower bound on the reheating temperature after modulus 
decay to thermalize the neutrino background for successful BBN (Tr > 1.2 MeV). Here we adopted 
= 3 and = 1. We also assumed that the abundance of the gravitinos produced through 
the reheating process after the primordial inflation are negligible. This corresponds to the case 
that the reheating temperature after primordial inflation is sufficiently low, = — 10^^ GeV, 
or J < Jcr- 
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0.1 1 10 102 

m3/3 (TeV) 

FIG. 2. Same as Fig. 1, but for various reheating temperatures after inflation, Tfj^inf = 10^^, 10^^ 
and 10^6 GeV, which are noted close to the corresponding lines. The extent of oblique lines coincides 
exactly with those excluded from the constraints by BBN and LSP in Fig. 1, which correspond 
to the limiting cases that Tr i^f is sufficiently low, T^^inf = — 10^^ GeV. The shadowed region 
corresponds to the excluded area for the highest possible reheating temperature, TR^inf = 10^^ 
GeV. The shadowed region corresponds to the excluded area for the highest possible reheating 
temperature, T^^inf = 10^^ GeV. Here we have assumed that = Mq but other cases can be 
easily read off from the fact that each line in the figure is determined by the value of J. 
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FIG. 3. Same as Fig. 1, but for Bh = I 
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TABLES 



"^3/2 




= 0-10^2 


TR,inf = 10^^ 


Ti?,inf 


= lO^'^ GeV 


0.1 TeV 


2 X 10^ 


- 3 X 10'^ TeV 


2 X 10^ - 


1 X 10"^ TeV 


2 X 10^ 


- 2 X 10^ TeV 


0.3 TeV 


8 X 10^ 


- 3 X 10'^ TeV 


8 X 10^ - 


1 X 10'^ TeV 


excluded 


1 TeV 


1 X 10^ 


- 1 X 10^ TeV 


1 X 10^ - 


8 X 10^ TeV 


1 X 10^ 


- 2 X 10*^ TeV 


3 TeV 


3 X 10^ 


- 2 X 10^ TcV 


3 X 10^ - 


1 X 10^ TeV 


3 X 10^ 


- 3 X 10^ TcV 


10 TeV 


2 X 10^ 


- 2 X 10^ TeV 


2 X 10^ - 


2 X 10^ TeV 


2 X 10^ 


- 7 X 10^ TeV 


30 TeV 


5 X 10^ 


- 2 X 10^° TeV 


5 X 10'' - 


2 X 10^° TeV 


5 X 10^ 


- 4 X 10^ TeV 


100 TeV 


2 X 10^ 


- 2 X 10^° TeV 


2 X 10^ - 


2 X 10^0 TeV 


2 X 10^ 


- 4 X 10^ TeV 



TABLE I. Allowed values of moduli mass m<^ for various 7713/2 and TR^inf for Bh = 10 



^3/2 


7R,inf= O-IO12 GeV 


7R,inf = 10^4 




TR,inf = 10^6 GeV 


0.1 TeV 


2 X lO'' - 2 X lO'^ TcV 


2 X lO'' - 1 X 10^ 


TeV 


2 X 10^ - 1 X 10^ 


TeV 


0.3 TeV 


1 X 10^ - 2 X 10*^ TeV 


excluded 




excluded 




1 TeV 


excluded 


excluded 




excluded 




3 TeV 


excluded 


excluded 




excluded 




10 TeV 


1 X 10^ - 3 X 10*^ TeV 


1 X 10^ - 3 X 10^ 


TeV 


1 X 10^ - 2 X 10^ 


TcV 


30 TeV 


5 X 10^ - 2 X 10^ TeV 


5 X 10^ - 2 X 10^ 


TeV 


5 X 10^ - 4 X 10^ 


TeV 


100 TeV 


2 X 10^ - 2 X 10^° TcV 


2 X 10^ - 2 X 10^° 


TeV 


2 X 10^ - 4 X 10^ 


TeV 



TABLE II. Same as Table 1, but for Bh = 1. 
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